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We compute the first order chiral phase transition for an instanton motivated quark model with a 
local six-quark interaction. In order to compare different solutions of the gap equation we compute 
the bosonic effective action - a two particle irreducible free energy functional. We find that the first 
order transition ends for a critical current quark mass, with continuous crossover for larger quark 
masses. Furthermore, we investigate different possible order parameters, including a color octet 
condensate. We also compare our formalism with mean field theory. 
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I. INTRODUCTION 

The understanding of strongly interacting fermion sys- 
tems is one of the great challenges for theoretical physics, 
ranging from condensed matter systems and ultracold 
atoms to strong interactions in particle physics. Non- 
perturbativc methods need to be explored in order to deal 
with strong effective interactions, collective phenomena, 
condensates in competing channels and different charac- 
teristic physics at "microscopic" and "macroscopic" dis- 
tances. In this note we deal with the bosonic effective 
action based on a two-particle irreducible formalism and 
the Schwinger-Dyson equations in a situation where the 
dominant interactions involve more than four fermions. 
We also compare with mean field theory. 

Indeed, for some physical systems the higher order 
fermion interactions play a crucial role. As a concrete 
example we discuss the instanton interaction in quantum 
chromodynamics (QCD) with three flavors of quarks. 
The violation of the axial J7(l)yi-symmctry by instanton 
effects induces an effective six-quark-interaction which 
involves all three flavors . It has been speculated that 
this instanton mediated interaction may become domi- 
nant at a characteristic momentum scale below IGeV. 
We explore here a model where the instanton mediated 
interaction dominates the effective theory at some given 
scale A. This scale will act as an ultraviolet cutoff for 
the fluctuations with momenta < A^ which induce 
condensates and spontaneous chiral symmetry breaking. 
We concentrate on the pointlike limit of this interaction 
and neglect all other effective interactions between the 
quarks. Wc do not believe that this model is a sufficient 
description for real QCD but it points to some character- 
istic features and serves as an interesting demonstration 
for the more formal part of this work. 

First investigations of strongly interacting fcrmionic 
systems are often based on Mean Field Theory (MFT) 
or lowest order Schwinger-Dyson equations (SDE) 



0, ^- For example, the recent studies of color su- 
perconductivity Ulllallllllllllll are mainly based 
on one of these methods. Both methods allow for a com- 
putation of the order parameter in systems which ex- 
hibit spontaneous symmetry breaking (SSB). However, 
while the SDE approach leads directly to the gap equa- 
tion, the MFT approach provides naturally a free en- 
ergy functional for the bosonic composite degrees of free- 
dom introduced by partial bosonization via a Hubbard- 
Stratonovich transformation [T^ ITU . The minima of the 
free energy are determined by the field equation which 
corresponds to a type of gap equation. Knowledge of the 
free energy functional becomes necessary if the gap (or 
field) equation allows solutions with different order pa- 
rameters and the free energy for the different solutions 
has to be compared. 

Within the SDE approach the reconstruction of the 
free energy functional from the gap equation is not triv- 
ial since the information about the minima (as expressed 
by the gap equation) may be insufficient in order to find 
the whole function. For example, it is not sure that 
the method used in for the case of color supercon- 
ductivity always works. From this point of view MFT 
seems superior to SDE since it directly provides a free 
energy. Unfortunately, MFT has also a severe disadvan- 
tage: partial bosonization is not unique and the results 
of the MFT calculation depend strongly on the choice 
of the mean field (see ^3 for the example of the Hub- 
bard model where this has drastic consequences). Since 
partial bosonization is an exact procedure physical re- 
sults should be independent of this choice. Indeed, the 
ambiguity is due to the approximation used in MFT 
(only fcrmionic fluctuations are taken into account) and 
is cured by more sophisticated approximations |l8l |. Be- 
ing a purely fcrmionic formulation SDE's do not suffer 
from such an ambiguity^. 
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^ For some reviews on SDE's and more sophisticated approxima- 
tions see e.g. |^|^|^. 

^ Moreover SDE's are one-loop exact even in the lowest approxi- 
mation, while for MFT this is not necessarily the case. This is 
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Hence, we want to find a functional which has the SDE 
as its equation of motion, and which can be interpreted as 
a free energy. For four fermion interactions such a func- 
tional can be related to the two particle irreducible (2PI) 
effective action JJ),^, 21_,.22J. The 2PI effective action is 
a functional of fields and propagators F'^^-'^^^ ["^i G] . How- 
ever, for a purely fermionic system, all the information 
is already contained in F^^^^' [0, G] where the "fermionic 
background field" is set to zero. Hence F^^^^^ [0, G] de- 
pends only on the bosonic variable G, and therefore we 
wiU call it the Bosonic Effective Action (BEA) |[23|. In a 
suitable version the BEA is a free energy functional and 
the condition for its local minima precisely corresponds 
to the gap equation. 

The work of (2^ has concentrated on effective four 
fermion interactions. In the present paper we want to 
deal with more general multi- fermion interactions. An 
example of such a higher order interaction is the six- 
fermion interaction generated by instantons in the case 
of three flavors and three colors Q, S HI El IH 113 ■ 
This interaction is U(1)a -anomalous and solves the fa- 
mous [/(l)A-problem j23| in QCD. In the simpler case 
of two flavors the instantons mediate a four fermion in- 
teraction which has already been studied extensively for 
chiral symmetry breaking [s^, |^ and was consid- 
ered in the early works on color superconductivity e.g. 
d 01 . This investigation is generalized here to three 
light flavors of quarks which are perhaps closer to realis- 
tic QCD. We also perform a more systematic discussion 
of the free energy functional. We observe that the effec- 
tive interaction generated by the instantons does not only 
lead to interactions between color singlet effective quark- 
antiquark degrees of freedom which arc finally associated 
to the usual spontaneous chiral symmetry breaking. It 
also produces effective interactions between color octets. 
In a model that goes beyond the six quark interactions 
the instanton effects may lead to the possibility of octet 
condensation and spontaneous "color symmetry break- 

ing" mmiii. 

The paper is organized as follows. In sects. IhI and Hill 
we introduce the 2PI formalism in the language of the 
Bosonic Effective Action (BEA). In sect. IIVI we make 
a short comparison to MFT and hint to some possible 
improvements. We explicitly calculate the BEA for the 
six- fermion instanton interaction in sect. Ivland use it to 
study chiral symmetry breaking. Finally, in sect IVII we 
investigate color octet condensates. Sect. IVIII summa- 
rizes our results and conclusions. 



effective action [23 and generalize it to multi-fermion in- 
teractions. In order to simplify the presentation we sum- 
marize all indices of the fermionic field in ^a- Here the 
index a contains all internal indices (spin, color, flavor 
etc.) as well as position or momentum. Furthermore it 
also differentiates between and ^p. In these conventions 
the partition function reads 

= J Vljj expirjatpa + ^japlpaiJlS " Sint[lp])- (1) 

All terms quadratic in tp are associated to bosonic sources 
jap and we investigate multifermion interactions (n even, 
n > 4) 



(2) 



The generating functional of the IPI Greens functions 
in presence of the bosonic sources j is defined by a Legen- 
dre transform with respect to the fermionic source term 
77: 



(3) 



where 



dW 



W = \nZ[T^,j], = (4) 

The fermionic effective action F^ can also be obtained 
by the implicite functional integral^: 



^f['>P,3] = -In J Vijje^p{ria-ipa - Sjl-ij) + Tp]), (5) 

This form is especially useful to derive the SDE. Taking 
a derivative with respect to ^ one finds 



- -Jl3a2Va2 + F Va2 



(6) 



where 

= (n-2)(n-4)---2 (7) 



II. BOSONIC EFFECTIVE ACTION 

In this section we briefly summarize the Schwinger 
Dyson equations 0, Q and the fermionic and bosonic 



(2) 

and Tp' denotes the second functional derivative. Here 
Za3...a„ summarizes all terms containing third and higher 
derivatives of F. (The corresponding diagrams involve 



directly related to the ambiguity. 



^ Note that in this formula tp is shifted such that (tp) = and 



STf 



depends on ip. 
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at least two vertices.) Taking another derivative with 
respect to ipa and evaluating at V' = we find the SDE: 



A 



(n) 

Q/3Q3...a„ 

'F'n 



(8) 



(2)n-1 



In this paper we are only interested in the lowest order. 
Therefore, from now on, we neglect Z, i.e. terms with 
more than one vertex. 

The "Bosonic Effective Action" (BEA) 22J, is defined 
by a different Legendre transform with respect to j for 
?7 = 0: 



TbIG] = -W^[0,j]+jG, 



^al3 — —. — F )c 



djaP 



/3 



(9) 

■jaf3. (10) 



Since Tf is an even functional of ip the BEA contains 
the same information as F^?. Indeed it is related to F^? 
by means of functional differential equations like Eq. 
(|10|l . Using this relation we can conveniently write the 
SDE as 



(n) 



Go,.^o,^ • • • Gan_ia.n • (11) 



Eq. H10() then yields a differential equation for Tb 



A 



in) 

q/3q3 .. .a„ 



^ n 

By integration"* one finally finds 



1 \V"-> 

Tb = -Tr In G + V ^i-^G„,„, • • • G„„_,«„ , (13) 
'> ^ nFn 



the BEA at "one- vertex order" . Actually it is sometimes 
convenient to introduce an auxiliary effective action 



^j[G,j] = Tb - -japGap 



(14) 



such that the physical propagator corresponds to the 
minimum of Tj (cf. Eq. 1)10(1 \ The functional Tj will 
play the role of a suitable free energy (see for details). 



III. BEA FOR LOCAL INTERACTIONS 

In the following we want to consider local interactions. 
For clarity we now write x (or momentum p) explic- 
itly and use latin letters for the remaining indices. The 



Note that in our notation 



— = Sa^/Sax — S„xSa 



standard procedure would be the insertion of the ansatz 
Gabi^^y) = -jabix,y) + Aab{x)6{x - y) into Eq. ||nil- 
This would yield the SDE for the local gap A. Since the 
BEA (O is related to the SDE ^ by differentiation 
with respect to G it is not clear, however, that an effec- 
tive action functional depending on A can be obtained by 
integration with respect to A. In presence of several pos- 
sible gaps this would require suitable "integrability con- 
ditions" for the system of gap equations. This difficulty 
can be avoided if we follow the construction presented in 
[2^ and start directly from the approximate BEA H13|l. 
With 



we have 



gab{x) = Gab{.X,x) 



iTrlnG+iTr(Gj) 



(15) 



(16) 



X! "np"" ^'^"-1"" 



For this relation it is essential that the interaction is 
strictly local. Furthermore, we can use the locality of 
the interaction in order to write Eq. Ulljl in the form of 
a local gap equation 

Gal^^^y) = -jab{x,y) + Aab{x)6{x-y). (17) 

We will evaluate the functional Tj [G] for Gap taking val- 
ues corresponding to Eq. I|17|l . This is actually a restric- 
tion to a subspace of all possible G. However, locality 
tells us that the extremum (solution of the SDE) is con- 
tained in this subspace. 

Using j = — G^* + A we find (up to a shift in the 
irrelevant constant and using Aab{x, y) — Aab{x)6{x—y)) 



iTrln(-j + A)-i 



Aab{x)gabix) 



"^p"" 9a^a2ix)■■■ga„_^aAx)■ (18) 



For the search of extrema of Tj it is actually convenient to 
treat A and g as independent variables. The extremum 
of Tj[g, A] then obeys 



gr,-[g,A] 
dA 



0, 



gr,[g,A] 
dg 



0. 



(19) 



Evaluating the derivative with respect to A we recover 
the inverse of Eq. H17|) for x = y, 



gabix) = i-j + A)J{x,x) = g[Aix)]. 



(20) 



Inserting this functional relation into Eq. Hll|) leads to a 
gap equation for A. In case of a six-fermion interaction 
this takes, however, the form of a two-loop equation. 

For n-fermion interactions with n > 4 it is more ap- 
propriate to go the other way around and first take a 



4 



derivative with respect to g. We obtain 



A 



Fn 

^ab[g{x)] 



(21) 



which is precisely the value of the gap in Eq. . Insert- 
ing A[g] into H18() we find the effective action depending 
on g 

r,[5] - -iTrln(-j + ALg])-i / A,,[g]ix)gabix) 



Searching for an extremum yields 



or 



J [91 



dg 



{i-j + A[g])-' - g} 



dA[g] 
dg 



= 0. 



(23) 



For ^ 7^ Eq. ^ indeed corresponds to the SDE 
(jni, i.e. 



gabix) = i-j + A[g])J{x,x). 



(24) 



Eq. H24I) will be our central gap equation. We should 
point out that possible cxtrcma of Tj [g] corresponding 
to ^ = are not solutions of the gap equation and 
should be discarded. Finally, we also have 



dg 



dr,[g,A[g]] 



dg 

dT,[g,A[g]] dT,[g, A[g]] dA[g] 

dg 

dT,[g,A[g]] dA[g] 



(25) 



dA 



dg 



dA 



dg 



Only as long as ^^^^ 7^ is fulfilled we can conclude that 
a solution of Eq. H23(l fulfills both extremum conditions 

The procedure proposed here is quite powerful if 
Trln(— j + A) can be explicitly evaluated as a functional 
of A. Then Tjlg] allows not only a search for the ex- 
tremum (discarding those with ^ = 0) but also a simple 
direct comparison of the relative free energy of different 
local extrema. This is crucial for the determination of 
the ground state in the case of several "competing gaps" . 

Furthermore, the formula H22(l for Tj (and the corre- 
sponding field or gap equation) is now a one loop expres- 
sion. This "one-loop" form of the equation of motion and 
the effective action is very close to what we would expect 
from MFT (cf. also the next section). In contrast to the 
standard SDE, which is only an equation of motion, we 
can use Eq. (|22|l to compare the values for the effective 
action at different solutions of the equation of motion 
(|23|l . providing us with information about the stability 
of a given state. 



The one vertex approximation (|21|l , (|22|l to the bosonic 
effective action is the central tool of this paper 



Vig{x)) 



V{g{x)) 



(2^ 



trln(-j+A[g]), (26) 



E(n — 2)Aii^..a„ . ^ / N 
7^ 5aia2 (a;) • ■ ■ 5a„-ia„ (.x) . 



(27) 

The first term may be called "classical part" and we have 
written the second "one loop term" as a momentum inte- 
gral with tr over internal indices and — A involving the 
Fourier transform of the gap (|21|l (see for details). 
As mentioned already we can use T j in order to find and 
compare the "local minima" of the free energy. Some care 
is needed, however, for the computation of susceptibilities 
or effective masses. Indeed, we have to be careful when 
considering Eq. (|22|l at points which are not solutions 
of Eq. Going step by step through the procedure 

above, we find that if we are not at a solution of (|23|) we 
do not necessarily fulfill the ansatz Ijl?!) . Therefore, at 
these points we are mathematically not allowed to insert 
the ansatz into Eq. \i'6\ . So, strictly speaking Eq. H22|) 
only gives the value of the effective action at the solu- 
tion of the equation of motion^. This is already much 
more than what we get from the standard SDE. Going 
beyond this we would also like to interpret Eq. (|22|l as 
a reasonable approximation in a small neighborhood of 
the solution to the equation of motion. Remembering 
g{x) = {ip{x)ip{x)) it is suggestive to interpret 5 as a 
bosonic field. Eq. H21|) gives the (non-linear) "Yukawa 
coupling" of g to the fermions, i.e. the relation between 
the gap and the bosonic field. Thus the term ^ Tr In is 
the contribution from the fermionic loop in a background 
field g. The second term in Eq. (|26|l can then be inter- 
preted as the cost in free energy for different background 
fields g. This interpretation allows us to use H22|l to cal- 
culate the mass and the couplings of the bosonic field g 
approximately. 



IV. COMPARISON WITH MFT 

We refer here to MFT as used in most computations 
and obtained by partial bosonization (for a more detailed 
description of the procedure and its ambiguities see, e.g. 



^ An alternative would be to choose the gap A as the "bosonic 
field". Inserting Eq. ^ into Eq. El 

we could calculate a 

functional r[A]. However, as one can check there are two draw- 
backs. First, even for four-fermion interactions, r[A] is usually 
unbounded from below when considering A 00. Second, in 
the case of a large four-fermion coupling the "stable" solution of 
the field equation is usually a local maximum. 
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[13 )• Using the identity 



1=JV X'0exp(TO(-0)*) =7V / V(l)exp{m{ip-{jj-(j>)^) 



(28) 

we can introduce bosonic fields corresponding to 
fermionic bilinears 4>ab(x) = {ipa{x)ipb{x)) . The coefB- 
cient m wiU be chosen such that the insertion of Eq. 
(|28|l into the functional integral iQJ cancels the purely 
fermionic interaction Si-at [ip] ■ We have not displayed the 
indices in Eq. H28|) . Usually one restricts the choice 
to one or several particular index pairs (a, b) (or linear 
combinations thereof). This corresponds to the freedom 
in the choice of the mean fields (j)ab{x). The partially 
bosonized form of Eq. (O becomes 



+1 (") 

'^aibi...a,i 

7 7 



•'^ n 

X \ 4>aibAx) ■ ■ ■ (t^ar^bn. i^) 

2 2 

-001 (a;)^/'&i {x)4>a2b2 [x] - ■■ (/)a„ 6„ {x) 

+ (1^2) + (1^3) + --- + (1 
+ --- + 0(?/;"-2) 



(29) 



and the functional integral for the partition function 
has now to be performed over bosonic as well as fermionic 
variables. We note that the coefficients m are only partly 
determined by A, i.e. 



/»ai ...an 



~\~ ...an 



(30) 



the "gap" A differ. Furthermore, we note that the coef- 
ficients m are not completely fixed due to the presence 
of the arbitrary symmetric part S in Eq. 13()ll . Results 
can depend on the choice of 5* (see e.g. |17k 18]). This 
is the so called "Fierz ambiguity" because the addition 
of S corresponds to a (generalized) Fierz transformation 
in the fermionic language. Of course, further considera- 
tions as e.g. the stability of the initial bosonic potential 
might reduce the freedom in the choice of S somewhat. 
But, as the example of 0| shows, such restrictions are 
sometimes not even strong enough to get qualitatively 
the same phase diagram for all reasonable choices of S. 

Eqs. ifTHl) . (j^ . do not suffer from such an am- 
biguity since in the derivation of the SDE Q the co- 
efficient becomes antisymmetrized and symmetric terms 
drop out. In 18] it was shown that the inclusion of ap- 
propriate diagrams for the bosonic fluctuations cures the 
Fierz ambiguity for four fermion interactions and leads to 
the SD-result. We believe that this holds also for higher 
fermion interactions. Nevertheless, the inclusion of the 
bosonic fluctuations needs a substantial effort. We there- 
fore propose H22|l as a natural replacement of Eq. H31() 
which amounts to an "optimal" choice of m for many pur- 
poses. Moreover, allowing not only for constant but also 
for spatially varying g we can calculate the wave function 
renormalizations and masses of the bosons. Again, the 
BE A gives unambiguous results. 

Finally, let us stress again the intuitive argument for 
the closeness of the two approaches and explicitly for 
the similarity of Eqs. (|22|) and (|31|l : we have g{x) = 
{i]j{x)^j{x)) which is exactly what one has in mind as a 
" mean field" . 



V. CHIRAL SYMMETRY BREAKING FROM A 
THREE-FLAVOR INSTANTON INTERACTION 



Here is a sum of terms which are symmetric in at least 
one pair of indices. The condition (|30|l ensures that the 
partially bosonized Lagrangian is equivalent to the origi- 
nal fermionic one. Nevertheless, due to the anticommut- 
ing nature of the fermionic variables S gives a vanishing 
contribution to the purely fermionic part of the action 
and is therefore not fixed. 

Neglecting the terms ©(V''*) and performing the func- 
tional integral over the fermions provides us with the 
MF-effective action: 

r^^[0] = jv[<^{x))-\TTH-3 + m), (31) 



ma^b, = (-1)* 



'a2&2 ■ ■ ■ </'aaf),i H ^ (1 ^ Tt) 

5 2 2 2 



This form is strikingly similar to Eq. H26|l . However, 
the coefficients in the "classical potential" as well as in 



In this section we want to use the method described 
above to study chiral symmetry breaking in an NJL-type 
model with a six-fermion interaction. We consider the 
QCD-instanton interaction with three colors and three 
flavors 0, mm mm in the pointlike Hmit. The three 
flavor instanton vertex can be written in the following 
convenient form l27l 



c 



£010203 £f)i62b3 



~{R ^ L) 



(32) 
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Here are the Gell-Mann matrices with color indices 
acting as generators of the SU {3)c color group. The 
brackets ( ) indicate contractions over color and spinor 
indices. Within QCD the coupling constant C should be 
be calculated in terms of the running gauge coupling. 
However, already the one loop approximation involves 
an IR divergent integral over the instanton size. There- 
fore, one needs to provide a physical cutoff mechanism. 
To avoid this difficulty we treat C as a free parameter in 
an effective theory with ultraviolet cutoff A. Inspection 
of (|32|l tells us that this interaction is C/(l)^-anomalous 
with a residual Z3-symmetry. This is important because 
we cannot restrict ourselves to real condensates from the 
start. 

In order to extract the coupling A*^^^ we have to anti- 
symmetrize over flavor indices (a = 1 ... 3), color indices 
{i ~ I. . .5) Weyl spinor indices (a = 1,2), chirality in- 
dices (x = 1,2 = L,R) and the indices distinguishing 
between ijj and (s = 1, 2): 



A. 



Y^^ai 0203 ^a4a^aQ f^a2CK5 ^a^aQ 



(33) 



\Z 1 

i2*5 ^13*6 J 



^ [ '^X 11*^X2 1*^X3 1*^X4 2*^X5 2 (5x6 2 
X [ Ss-^2Ss22Ss32SsilSs^lSssi 

+ (^S 1 1 <5s2 1 <5s3 1 <5s4 2 <5s5 2 <5s6 2] • 

Here P denotes the sum over all 6! permutations of the 
multiindices rrij = {aj,ij,aj,Xj^Sj), j = 1...6, with 
minus signs appropriate for total antisymmetrization. 

As a first example we consider a flavor singlet, 
color singlet scalar chiral bilinear {a = ^''Pl'ipii, 



(j{Sxl^T2Ss2Stl - Sx2SrlSslSt2) 



-a*{6^2STlSs2Stl - SxlSr2SslSt2) 



(34) 



We evaluate 



^mm2m3nm5mQ 



(35) 



-6 



jm^m^gmsTae 



g- 



m^megm^m^ 



where we have used the fact that A is symmetric under 
permutations of the three ipij: bilinears. Exploiting the 



flavor, spin and color structure for H35() expresses the gap 
in terms of a 



(T'^{5xl5T25s25tl - 5x2^Tl^sl5t2) 



-(T*^((5x2'5rl(5s2'5tl - 6xi5t25si5t2) 



(36) 



For the classical potential one obtains 
20 



V{a) 



(37) 



For a background a which is constant in space we obtain 
the effective potential U {x — q^) 



U{a) 



dx x[lTi{x+\Mg\^)] + V{a) 
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(38) 



Here rriq is a current quark mass which we take to be 
equal for all quarks. The integral is, of course, di- 
vergent. Our UV regularization simply cuts it off with 
X < A^. Measuring all quantities in units of A we can 
put A = 1. 



A. The chiral limit — 

Let us now search for cxtrema of U{a)^. Inspection 
of U{a) tells us that it is invariant under the combined 
operation C— s-— C, cr— >— cr. This allows us to restrict 
our analysis to positive ^. In the chiral limit it is useful 
to parametrize 



a = \a\ exp(zQ;). 



From this one finds 



40 



C^(kU) = yCk|'cos(3a) + /(|(T|), 



(39) 



(40) 



where / is a function determined by the integral in Eq. 
(|38|l . We can see that the only a-dependence comes from 
cos(3q;) which is the explicit manifestation of the Z^- 
symmetry. It is clear that extrema can only occur at 
a — n G Z. Using the Zs-symmetry we can restrict 
ourselves to a = 0, tt or restrict ourselves simply to real 
a. 

Taking all this into account we find up to three so- 
lutions (cf. Fig. By symmetry cr = is a solu- 
tion for all values of the coupling. Going to larger cou- 
plings we encounter a point Cs where we have two solu- 
tions. For even larger couplings there are three solutions 



Since ' =^ for all cr 7^ we do not need to worry for non- 
trivial solutions to be spurious. In addition cr = is always a 
solution in the chiral limit. 
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FIG. 1: Effective potential for various values of the coupling 
constant increasing from the topmost line C, — 3000 to the 
lowest line ( = 4200. The second line (long dashed), for 
<^s ~ 3350, corresponds to the critical coupling for the onset 
of non-vanishing solutions. The third line is for ( = 3600 
while the fourth (short dashed) denotes the onset of SSB at 
CssB ~ 3900. The horizontal line indicates the value of U{a) 
at the trivial solution a = 0. 



= o-Q < 0-1 < 0-2- We know that U{ai) > U{ao = 0) 
therefore ai is not the stable solution. As can be seen 
from Figs. QElthere is a range Cs < C < Cssb where there 
exist non-trivial solutions to the SDE but no spontaneous 
chiral symmetry breaking occurs because U{a2) > U{Q). 
The state with cr ^ is mctastable. This holds as long as 
Cs < C < Cssb, with Cs the point of "spinodal decomposi- 
tion" . A first order phase transition towards a state with 
nonvanishing chiral condensate a occurs as C, is increased 
beyond Cssb- Here the critical value Cssb for the onset of 
spontaneous chiral symmetry breaking denotes the value 
of C for which the two local minima become degenerate, 
U{<J2) — f/(0). We point out that in order to calculate 
Cssb we need to know the value of [/, i.e. information 
beyond the SDE. 

In Fig. 121 we have plotted the mass gap versus the 
six-fermion coupling strength. We observe the first order 
phase transition. This may be expected generically due 
to the "cubic term" in the classical potential (|S7|l . 

Finally, we would like to remark that in general the 
free energy does not follow from a naive integration of 
the SDE with respect to the gap A. We may call the re- 
sult of a naive integration of the gap equation the "pseudo 
potential" C/(A), or, more generally, the "pseudo-free en- 
ergy" r[A]. In general r[A] is not equal to the the BE A 
rj [G], not even at solutions of the SDE. Indeed, as can 
be seen from Fig. [S] the results for physical quantities 
like the effective fermion mass can differ. The underlying 
reason for this is that the gap A is not the correct inte- 
gration variable. The SDE is obtained by a G-derivative 
of the BEA functional rj[G]. Therefore, in order to re- 
construct rj[G], we have to integrate with respect to G. 
As can be seen from Eq. (|21|l A is, in general, not even a 
linear function of G. Simple integration with respect to 



FIG. 2: Constituent quark mass \Mq\ (gap) versus the 
strength of the six-fermion interaction C, in the chiral limit 
niq — 0. The thick line corresponds to the solution with small- 
est action. The dashed line reflects the largest non-trivial so- 
lution (T = (72 . Finally, the thin line is obtained by minimizing 
a "pseudo effective potential" U (A) obtained by direct inte- 
gration of the SDE with respect to the gap A. We find three 
special couplings, (^s for the onset of non-trivial solutions to 
the SDE, (^ssB for the onset of SSB, i.e. a non-trivial solution 
that has lower action than the trivial solution and Cssb where 
the lowest extremum of U{A) becomes non-trivial. In our ap- 
proximation we obtain a first order phase transition. We also 
see that the "true" phase transition at Cssb differs from the 
value obtained by a naive integration of the SDE. 

A therefore neglects the Jacobi matrix, which is a non- 
trivial function of A for interactions more complicated 
than a four-fcrmion interaction. 



B. Non- vanishing current quark masses 7^ 

The non-vanishing current quark mass explicitly 
breaks the residual Za-symmetry. The effective poten- 
tial C/(|cr|, a) does no longer depend on cos(3q) only, and 
we have to look at the complete complex cr-plane for pos- 
sible extrema. Moreover, for ruq = 0, U{(7) is completely 
symmetric under C —* ~Ci'^ ~* ~cr- Therefore, we could 
restrict ourselves to C > 0. For 7^ we need to add 
the transformation — > — m^. We can still restrict our- 
selves to positive C but we need to consider both positive 
and negative m^. 

In the case of rUq 7^ we still encounter an extremum of 
U{a) at (7 = 0. However, in this case it is not a solution of 
the SDE. It is a spurious solution due to ^^jr^ = 0. The 
difference to the chiral limit is that the derivative of the 
effective potential U{a) now has only a simple zero while 
in the chiral limit it has a twofold zero. After dividing 
the field equation by ^^^5^ a simple zero remains only in 
the chiral limit, giving a solution of the SDE. 

Although chiral symmetry is now broken explicitly we 
can still observe a first-order phase transition signaled 
by a jump in the fermion mass. The critical coupling for 



8 




1. 

0.8 

0.6 

0.4 

0.2 



0.02 



0.04 



0.06 



0.08 



FIG. 3: Dependence of the critical coupling ("ssb on the cur- 
rent quark mass (solid line) and jump in the fermion mass 
at the phase transition |5Mq| (dashed line, scale on the right 
hand side). We observe that there exists a critical ~ 0.076 
above which the phase transition is replaced by a continuous 
crossover. 



the phase transition depends on ruq as depicted in Fig. 
13 The critical hne ends at niq = mq^crit ~ 0.076. For 
niq > TOq crit wc havB no first order phase transition in 
our approximation but rather find a continuous crossover. 
We have also plotted in Fig. Othe jump in the constituent 
quark mass \SMg\ between the phases with low and large 
\a\ at the critical coupling Cssb- It vanishes at the end 
of the critical line at Wg^crit- 



VI. COLOR-OCTET CONDENSATION 

In the preceeding section we have considered only one 
direction in the space of all possible g resulting in a phase 
diagram for chiral symmetry breaking. Let us now con- 
sider the more general case where we also allow for a 
non-vanishing expectation value in the color-octet chan- 
nel, more explicitly in the color-flavor locked direction 



:SabSijdap 
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A condensate in this direction would yield a very inter- 
esting phenomenology for QCD where the gluons acquire 
a mass via the Hi ggs mechanism and are associated with 
the p- mesons |32l l33j . A pointlike six-quark interaction 
seems too simple in this case for a realistic model of QCD 
- it does not account for the important infrared cutoff of 



the instanton size due to the effective gluon mass 
Nevertheless, we find a study of a pointlike interaction 
an interesting preparation for the treatment of a more 
realistic model. 

Following the outline of the previous section we obtain 
the gap 
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and the classical potential 



v{<j,x) = C 



20 



;(^X 



-,X'' + c.c. 



(42) 



(43) 



9" 18" 54^/6 

This results in the effective potential 

C/(a,x) = V[a,x] (44) 

- ^ I dx x[S\n{x + iMsH + hi(x + jAfip)] 
J57r^ Jx 

where the singlet and octet masses for the "constituent 
quarks" associated with the low mass baryons [s^ [s^ 
are given by 

Ml = lc(12a2-\/6ax-x')+"i9 (45) 

= J-C(192a2 + 2N/6ax-7x')+m,. (46) 
864 

We start by a search of extrema of U for real a and 
X- In the chiral limit every point on the line x = — 2V6cr 
(Ml = and Mg — 0) has the same value of J7 = 0, and 
both derivatives with respect to a and x vanish. How- 
ever, this is one of the spurious solutions mentioned in 
Sect, mil where ^ vanishes. (In our case A = (Mi, Afg) 

and g — (cr, x) are two component vectors and ^ stands 
for the Jacobian.) Direct insertion into Eq. (|24|l (A (5) = 
0) shows that on this hne only the point (cr, x) — (0, 0) is 
a true solution to the SDK. Restricting both a and x to 
be real we have not found a solution of the gap equation 
with X ^ 0. Thus, we have not identified a solution for 
which the condensate would break color symmetry but 
not parity. 

For the most general case of complex a and x things are 
considerably more difficult since we now have to search 
for an extremum of a potential which depends on four real 
parameters. We checked several values of the coupling 
constant. So far we have not found a solution which has 
a lower free energy than the minimum of the free energy 
for vanishing octet condensate x = 0. 

Still, we would like to point out that the potential is 
unbounded from below in various directions, including 
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those with X 7^ 0- Therefore, a physical cutoff mechanism 
hke the one discussed in [23 or a different approxima- 
tion to the "classical action" which makes the potential 
bounded from below may provide additional solutions. In 
this context we stress that the instanton interaction l|32() 
should not be confounded with realistic QCD. The gluon 
fluctuations have been omitted here. Perhaps even more 
important, the effective gluon mass for nonvanishing x 
should lead to an effective x-dependence of the coupling 
C HI]. Also the flat direction"^ of t/ for x = -l^a re- 
mains intriguing. So far we have not yet understood the 
reason why U , Mi and Mg all vanish simultaneously on 
this "line" . A small additional contribution to U (for ex- 
ample from the gauge boson fluctuations) could lift this 
degeneracy and lead to a minimum of U somewhat away 
from this line, such that ^ does not vanish anymore. 

Finally, let us compare the result of our SD calculation 
with the MFT result analogous to the computation in 
1351. However, we use here the corrected instanton vertex 
of |27| (without the cutoff mechanism for the instanton 
interaction considered there) and no other interaction. 
We apply the formulae of sect. llVl even though in our case 
the integral (|28|l is not well defined. Adopting the same 
normalization for the MP as for the condensates g in the 
SD calculation and taking the mean field as suggested by 
the brackets in (|32|l without further Fierz transformation 
the "classical MP potential" is 

^(CTmf, Xmf) = -C(o-MF + ^CTmfXmf) + C-C. (47) 

Por the "one-loop" Potential Umf we recover the form 
H44|l but with masses 

Mr = -C(4f + J/I'^mfXmf + ^xL) (48) 

This quite different from our SD result Eqs. (g3J). In 

particular the x'^-term in the classical potential is absent 
in the MFT calculation. Moreover, the sign between the 
(T'^-term and the trx^ is different, too. This demonstrates 
that the difference is more than an overall normalization 
of the potential or the fields. On the one side this high- 
lights once more the importance of the Fierz ambiguity. 
On the other side it is not obvious if a suitable mean 
field formulation exists at all which would reproduce the 
results of the SDK. In view of all the problems of MFT 
we would like to argue that the SDE is clearly superior 
for our problem. 



VII. SUMMARY AND CONCLUSIONS 

Integrating the lowest order Schwinger-Dyson equa- 
tion (SDE) for a multifermion-interaction we obtain the 
bosonic effective action at "one- vertex" -level. Within 
this approximation we find an " one-loop" expression for 
the SDE even in case of interactions involving more than 
four fermions. Although this gap equation is formally 
very similar to mean field theory, it does not suffer from 
the ambiguities of the latter arising from the bosonization 
procedure. We also propose a simple one loop formula 
for the free energy functional at the extrema which can 
be used in order to compare different local extrema. 

We apply our method to a six-fermion interaction re- 
sembling the instanton induced quark vertex for three 
colors and flavors. We compute the solutions of the gap 
equation and the minimum of the free energy in depen- 
dence of the coupling strength C,. For small current quark 
masses ruq the gap equation has several solutions, corre- 
sponding to different local extrema of the free energy. 
The free energy of the different extrema is compared by 
use of the one loop formula for the bosonic effective ac- 
tion. We find a first order transition to a phase with chi- 
ral symmetry breaking as the coupling of the six-quark 
vertex C increases beyond a critical value Cssb- Thereby 
the value of Cssb is larger than the value Cs , the minimal 
coupling for which solutions with a non-vanishing chiral 
order parameter exist. (Note that Cssb equals Cs only 
in the case of the perhaps more familiar second order 
transition.) The critical line in the space of the current 
quark mass niq and C ends at a mass mg_crit ~ 0.076 (in 
units of the UV cutoff near IGeV). For > mq^crit the 
phase diagram is characterized by a continuous crossover. 
We also investigate possible color octet condensates in 
the color-flavor-locked direction. In the approximation 
of a pointlike instanton induced six-quark interaction no 
phase with non-vanishing color octet condensate is vis- 
ible. In this respect our work should be considered as 
a starting point for a more realistic instanton induced 
interaction, where the dependence of the instanton solu- 
tion on the value of the octet condensate is taken into 
account. 

We conclude that the bosonic effective action pro- 
vides a practical tool for the understanding of fermionic 
systems where interactions involving more than four 
fermions play an important role. Beyond the gap equa- 
tion it provides for a free energy. In contrast to mean 
field theory the lowest order is well defined and gives an 
unambiguous answer. 
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This flat direction is also present for niq ^ 0. 
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